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The contribution of the QCD Pomeron to the process AA AA J/^ J/^ is discussed. We focus 
on the photon-photon colhsion, with the quasi-real photon coming from the Weizsacker - WiUiams 
spectrum of the nuclei. We calculate the cross section for this process considering the solution of the 
LLA BFKL equation at zero momentum transfer using a small t approximation for the differential 
cross section of the subprocess. Furthermore, the impact of non-leading corrections to the BFKL 
CO ' equation is also analyzed. In both cases the cross section is found to increase with the energy, 

' predicting considerable values for the LHC energies. Moreover, we compare our results with the 

, Born two-gluon approximation, which is energy independent at the photon level. Our results indicate 

(N , 

that the experimental analyzes of this process can be useful to discriminate the QCD dynamics at 
P3 ' high energies. 
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J I Understanding the behavior of high energy hadron reactions from a fundamental perspective within QCD is an 

important goal of particle physics. The behavior of scattering in the limit of high energy and fixed momentum 
Cn| I transfer is described in QCD, at least in situations in which perturbation applies, by the Balitskii-Fadin-Kuraev- 
^-H Lipatov (BFKL) Pomeron [1]. Attempts to test experimentally this sector of QCD have started in last years, mainly 
considering processes where specific conditions are satisfied, which minimize the contributions of the other mechanisms 
, competing with the QCD Pomeron and guarantee the validity of the perturbative QCD methods. In particular, the 
f~i ' virtualities of the gluons along the ladder should be large enough to assure the applicability of the perturbative 
Oh] expansion. The hard scale may be provided either by the coupling of the ladder to scattering particles, which contain 
a hard scale themselves, or by a large momentum transfer carried by the gluons. Furthermore, to distinguish the BFKL 
I from DGLAP evolution effects it is convenient to focus on processes in which the scales on both ends of the ladder 
I are of comparable size. Some examples are the measurements of forward jets in deeply inelastic events at low values 
of the Bjorken variable x in lepton-hadron scattering, jet production at large rapidity separations in hadron-hadron 
collisions and off-shell photon scattering at high energy in e+ e~ colliders, where the photons are produced from the 
leptons beams by bremsstrahlung (For a recent review of BFKL searches, see e.g. Ref. [2]). This last process presents 
some theoretical advantages as a probe of QCD Pomeron dynamics compared to the other ones because it does not 
involve a nonperturbative target [3-5]. Moreover, such reaction presents analogies with the process of scattering of 
two quarkonia, which has been proposed as a gedanken experiment to investigate the high energy regime in QCD [6]. 
In such a case, nonperturbative effects are suppressed by the smallness of the quarkonium radius. 

Other possibility for the study of the QCD Pomeron is the vector meson pairs production in 77 collisions [7] . At 
very high energies s ^ ~t, diffractive processes such as 77 neutral vector (or pseudoscalar) meson pairs with 
virtual or real photons can test the QCD Pomeron (Odderon) in a detailed way utilizing the simplest possible initial 
state. As in the case of the large angle exclusive 77 processes, the scattering amplitude is computed by convoluting 
the hard scattering pQCD amplitude for 77 — + qqqq with the vector meson wave functions. For heavy vector mesons, 
this cross section can be calculated using the perturbative QCD methods. First calculations considering the Born 
two-gluon approximation have been done in Refs. [7]. More recently, the double J/vP production in 77 collisions has 
been proposed as a probe of the hard QCD Pomeron [8]. There, the hard QCD pomeron is presumably the dominant 
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mechanism. Theoretical estimates of the cross sections presented in [8] have demonstrated that measurement of this 
reaction at a Photon Colhder should be feasible. 

In this paper we study the possibilities for investigating QCD Pomeron effects in a different context, namely in 
photon-photon scattering at ultraperipheral heavy ion collisions. In this case, the cross sections are enhanced since 
the 77 luminosity increases as Z^, where Z is the atomic number [9,10]. Here, we will focus our analyzes on double 
diffractive J/^ production in 77 collisions, with the photons coming from the Weizsacker - Williams spectrum of the 
nuclei. This process is unique since in principle it allows to test the QCD Pomeron for arbitrary momentum transfers, 
where the hard scale is already provided by the relative large mass of the charm quark. Here, our goal is twofold: 
to analyze the potentiality of this process to constraint the QCD dynamics at high values of energy and to provide 
reliable estimates for the cross sections concerning that reaction. We calculate the cross section considering suitable 
approximations for the 77 — > J/^* J/^ subprocess. In particular, we estimate it considering the two-gluon exchange 
(Born level), the leading order solution of the BFKL equation, as well as estimate the corrections on the energy 
dependence associated to the next-to-leading effects through the Pomeron intercept (for a review on NLO BFKL 
corrections, see e.g. Ref. [11] and references therein). Shortly, we have found that the two-gluon exchange mechanism 
gives rise to a constant 77 cross section at large c.m.s. two-photon energy W . To higher orders in perturbation theory, 
the iteration of gluons in the f-channel promotes this constant to logarithms, and the perturbative expansion of the 
cross section at high energy has the form, 



where Sq is the scale of the order of the charm mass squared, the summation represents the series of the leading 
logarithms on energy to all orders in the strong coupling constant a^, and the ellipsis stands for nonleading terms. 
The sum of the leading logarithmic terms is made by the BFKL equation at leading order, which predicts that 
cross section grows exponentially, a oc W'^"'^ , where ojp = 4 ln2 Nc as/ tt in the LLA solution. It has recently been 
demonstrated that the NLO corrections to the BFKL equation are large. The main effect is a reduced value of the 
so-called Lipatov exponent ivp. Since then the effects of higher orders have been studied for measurable processes, 
but the conclusions are not unambiguous [2] . This situation should be improved in the future with the next generation 
of linear colliders, such as the Japan Linear Collider (JLC) at KEK, TcV Energy Superconducting Linear Accelerator 
(TESLA) at DESY and CERN Linear Collider (CLIC), etc. However, until these colliders become reality is important 
to consider alternative searches in the current accelerators which allow us to constrain the QCD dynamics. 

This paper is organized as follows. In next section (Section II) we present a brief review of the ultraperipheral heavy 
ion collisions and the main formulae to describe the photon-photon process in these reactions. In Section III we consider 
the double J/^* production in 77 collisions and estimate this cross section considering distinct approximations to the 
QCD Pomeron dynamics. There, we calculate the total cross sections for the subprocess for each approximation and 
present our results for double J/^ production in ultraperipheral heavy ion collisions. Finally, we present a summary 
of our main conclusions in the Section IV. 



In ultraperipheral relativistic heavy-ion collisions the ions do not interact directly with each other and move essen- 
tially undisturbed along the beam direction [9] . The only possible interaction is due to the long range electromagnetic 
interaction and diffractive processes. Due to the coherent action of all the protons in the nucleus, the electromagnetic 
field is very strong and the resulting flux of equivalent photons is large. A photon stemming from the electromagnetic 
field of one of the two colliding nuclei can interact with one photon of the other nucleus (two-photon process) or can 
penetrate into the other nucleus and interact with its hadrons (photon- nucleus process). For a recent review, see e.g. 
Ref. [10]. In particular, the photoproduction of heavy quarks as a probe of the high density effects have been recently 
emphasized in Refs. [12], where the color glass condensate formalism [13] was used to estimate the cross section and 
transverse momentum spectrum. Similarly, the elastic photoproduction of vector mesons in ultraperipheral heavy 
ion collisions was studied recently in Ref. [14], demonstrating that this process can be used to constraint the nuclear 
gluon distribution, which determines the dynamics at high energies. For other studies of QCD dynamics in one-photon 
processes see Ref. [15]. Here, we will restrict our analyzes for the two-photon process and its potentiality to investigate 
the QCD dynamics. 

Relativistic heavy-ion collisions are a potentially prolific source of 77 collisions at high energy colliders. The 
advantage of using heavy ions is that the cross sections varies as Z^a'^ rather just as a^. Moreover, the maximum 
77 collision energy W^j is 27/i?^, about 6 GeV at RHIC and 200 GeV at LHC, where Ra is the nuclear radius and 
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7 is the center-of-mass system Lorcntz factor of each ion. In particular, the LHC wiU have a significant energy and 
luminosity reach beyond LEP2, and could be a bridge to 77 collisions at a future e+e^ linear collider. For two-photon 
collisions, the cross section for the reaction AA — > AAJ/'^ J/'^, represented in Fig. 1, will be given by 

(JAA^AAJ/<i^J/<5! = [ — niiui) [ — 712(^2) Cr^j^j/^j/^{W = Viw^) , (2) 
J 1^1 J (^2 

where the photon energy distribution n{cu) is calculated within the equivalent photon or Fermi- Weizsacker- Williams 
(FWW) approximation [16]. In general, the total cross section AA AAjj AAX, where X is the system 

produced within the rapidity gap, factorizes into the photon-photon luminosity and the cross section of the 77 
interaction, 

ctaa^aaj/^ j/<if{s) = j dr ^-fj^j/f jf-fis), (3) 

where r = s/s, s = W"^ is the square of the center of mass (c.m.s.) system energy of the two photons and s of the 
ion-ion system. The 77 luminosity is given by the convolution of the photon fluxes from two ultrarelativistic nuclei: 

dC^j (r) dx 



dT 



r dr 

= / -f{x)f{r/x), (4) 



where the photon distribution function /(.t) is related to the equivalent photon number n{uj) via f{x) — (E/lu) n{xE), 
with X = u)/E and E is the total energy of the initial particle in a given reference frame. The remaining quantity 
to be determined in order to proceed is the quantity f{x), which has been investigated by several groups (for more 
details, see e.g. [10]). Here, we consider the photon distribution of Ref. [17], providing a photon distribution which 
is not factorizable. The authors of [17] produced practical parametric expressions for the differential luminosity by 
adjusting the theoretical results, which reads as, 

dC^^_fZ^aYm ^{z), 0.05 < 2; < 5 



dr ~\ TT J 3t \ [In(if^)]^ z<Om 

where z = mR/^ = 2MR^, m is the mass of the system produced in the two photon collision, M is the nucleus 
mass, R its radius and ^{z) is given by, 

e(2) = ^A,e-^-, (6) 

i=l 

from an adjust to the numerical integration of the photon distribution, with an accuracy of ^ 2% in the referred region 
on z. The adjustable parameters are the following: Ai = 1.909, A2 = 12.35, As = 46.28, h = 2.566, 62 = 4.948, and 
63 = 15.21. 

The approach given above excludes possible final state interactions of the produced particles with the colliding 

particles, allowing reliable calculations of ultrapcriphcral heavy ion collisions. Therefore, to estimate the double J/^" 
production it is only necessary to consider a suitable QCD model for the double heavy meson production. In lines 
of the analysis presented here, it is worth mentioning that estimates for the double pion (light mesons) production 
in ultraperipheral heavy ion collisions have been presented in Ref. [18]. In addition, there it was found a negligible 
contribution of pomeron-pomeron (considered in the Regge approach) interaction for very heavy ions, whereas this is 
non- negligible for lightest ions. 



III. DOUBLE J/* PRODUCTION 

The process 77 J/'^J/'^ is a clean reaction testing the BFKL Pomeron physics and provides estimates for the 
heavy mesons production in photon induced processes. The hard scale involved in the reaction, that is the charm 
mass, and the presence of the photon as the initial state turns out it suitable for perturbative treatment. The non- 
perturbative content is provided only by the J/^ light-cone wave function, which is well constrained through the 
experimental measiucmcnt of its Icptonic width Fj/^,^ ;+;- . In the following we use the high energy factorization and 
the BFKL dynamics in order to perform estimates for the referred reaction. Our general formulae for the differential 
cross section and imaginary part of the scattering amplitude are given as follows [8] , 
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da{^^^J/^J/^)_\A{W\t)\^ 



dt 16 n 

-121, (ft ,._^Jl. „\ g, ,,_,_/'T/f/2 



Xm.4(l^,.) = y— (fe^^^,). , (8) 

where W is the center of mass energy of the two photon system and the photon-meson impact factor is denoted by 
At the Born level $ = $ and the reaction is described by the two-gluon (a bare approximation for the QCD 
Pomeron), which have transverse momenta q/2±k and where the momentum transfer is t = —q^. When considering 
the complete gluon ladder contribution, the quantity $ contains the impact factor and the gluon emission on the 
ladder, which is driven by the QCD dynamics. At the LLA level, the BFKL ladder contribution for the t-channel 
exchange provides the following expression for it, 

^^jMW\k,q) = J d'k' ^F{W\k,k',q)^^j/^{k',q) , (9) 



where 

^'yj/^{k,q) =C 



q^ ml + k^ 



(10) 



and FiW"^ ,k,k' ,q) is the solution of the LLA BFKL equation at f ^ 0. The Eq. (10) defines the impact factor 
in the nonrelativistic approximation. We have considered the following parameters for the further calculations: 
rric = = 1.55 GcV, C = ^a.em ctsip"^) ec | ttTOc ./,//*: with e,. = 2/3 and fj/q, = 0.38 GcV. 

In order to compute the total cross section for the process, we would need consider the LLA BFKL non-forward 
solution [19] (for a derivation in the momentum space, see e.g. [20]) in Eq. (9), put all together to obtain the 
amplitude and then integrate Eq. (7) over < |t| < oo. Further refinements can be considered. For instance, the 
infrared contributions by modifying the gluon propagator can be still investigated, as done in Ref. [21] for the Born 
case and as in Ref. [8] for the full resummation. Moreover, non-leading corrections could be also implemented in a 
numerical calculation, for instance in lines of Ref. [8]. Here, we follow a different procedure, since we are interested 
mostly in the energy dependence of the process. We will consider a small t approximation, providing accuracy enough 
for our purposes in the present work, in a such way that the total cross section can be written as. 



/ ./T 7/T^ 1 da{jj^ J/^J/^) 

otot (77 ^ J/-^) = -5 , , — 

-o J/* J/* «i 

with 

da{'y^^ J/^ J/^) \A{W'^,t = Q)\'^ 



(11) 



dt 16 TT 



exp {-Bji^ji^, ■ \t\) , (12) 



where B j/^ j/^, is the corresponding slope parameter. This quantity is not available from experimental measurements, 
instead of the slope Bjj^^ obtained from J/^E" photoproduction [22]. In the following we should estimate this quantity 
using information about the Born two-gluon approximation, where F = S^(k — k') . In this particular case the 
amplitude reads as, 

A rT/l/2 ,^ f d'^k %j/^{k,q) ^^j,^{k,q) 

which can be computed analytically if one considers a strong coupling constant not depending on k'^ . Using the impact 
factor given by Eq. (10) and performing the integration over gluon the transverse momentum, one obtains 

16C^ 



Now, we are able to use Eq. (14) in order to estimate the B parameter for the double J/^" production. This quantity 
can be computed from the scattering amplitude above by the small t approximation. 



which gives a slope parameter with vahie Bj/^, jj^ = ra~'^ . Moreover, we expect that this estimate has httle sensitivity 
on energy, since it is known that the hard Pomeron gives a low a'p (almost no shrinkage) [23]. That is, the phe- 
nomenology on J /'^ photoproduction has provided a'p ~ 0.1 [22], instead of the usual a'p = 0.25 from the hadronic 
analysis, where the corresponding slope is given hy B = bo + 2 a'p ln(s/so) in the Regge analysis. We have verified 
also the accuracy of our approximation comparing the analytical calculation, Eq. (15), with the small-f estimate given 
by Eq. (12). The comparison is shown in Fig. 2, for a fixed coupling constant asim^^ij,) ~ 0.289, in agreement with 
the estimates at Ref. [8]. It is worth mentioning that the typical values measured in elastic J/^* photoproduction stay 
in the range \t\ < 1.5 GeV^^, where our approximation produces the same result as the analytical case. 

Having the slope parameter B, and taking Eq. (11), one can estimate the total cross section considering the Born 
two-gluon exchange. It is simple to verify that the normalization for atot depends significantly on the strong coupling 
as- We have obtained that the cross section takes values of order a few pb's, having a lower bound if the typical 
HERA value ag ~ 0.2 is used as well as an upper bound for Q!s(m^). These estimates are corroborated by previous 
calculations [7], and some results are presented in the Table I, considering different values of as- The results are also 
consistent with Ref. [8], where the calculations were done beyond the leading log approximation. For instance, we 
have found aioi = 3.6 pb using the nmning as{k'^ + m^) and sq = 0.l6GeV^ (this is the parameter accounting for 
the infrared cutoff). The refined study of [8] gives a close result, atot ~ 2 — 2.6 pb. Our estimates are quite rehable 
in view of the approximations considered here. It should be stressed that other additional sources of uncertainty are 
the charm mass and the nonrelativistic approximation for the impact factor. 

In order to perform a LLA BFKL calculation, the following solution for the evolution equation in the forward case 
was considered (see the pedagogical reviews on Refs. [20,23]), 

1 1 / U^2 \ '^JP 

F(W^,k, k') = „ ., , — exp 

^ ^ ^/2n^ak'k" VHWVS) V W 

where 

ujp = ^Aln2, a=^28C(3) (17) 

and 1 + LOjp is the Pomeron intercept in the leading logarithmic approximation, which depends on a^. We have taken 
s = 1 GeV^. The results are sensitive to the choice for the intercept, providing an enhancement of the total cross 
section by one/two orders of magnitude in relation to the Born level at the considered energy range. For the further 
studies we have selected the results using as(/x^) = 0.22, typical at the HERA kinematical region. 

The total cross section is shown in Fig. 3, represented by the solid line, producing an effective behavior given by 
W"^^ and where A = 0.54. The range on W is the region possibly to be available at LEP2 and two-photon process 
in ultraperipheral heavy ions collisions at LHC. The behavior is considerably steep on energy and it is timely ask by 
the NLO corrections to the LLA BFKL equation or unitarity corrections to the LO calculation. Concerning the NLO 
effects, the correction looks like Wp^^ = top (1 — Nag), where = ^as/i^. The value A'' = 6.5, obtained in the 
literature [24] , turns out the corrections too large and even producing negative values if large Qg values are considered. 
The complete resummation of the non-leading effects, mostly collinear-enhanced contributions, is currently argued in 
order to provide reliable calculations (see Ref. [11] for a review). 

For practical purpose in our further investigation below, we have considered the modification of the Pomeron 
intercept in order to be consistent with the stable higher order resummation results [25], which give ujp ~ 0.3 at 
the HERA kinematic regime. Even a lower effective exponent is found concerning the NLO BFKL using the BLM 
scheme for the renormalization scale setting [26]. For a full calculation we would need the NLO impact factor, which 
is not completely known [27]. However, one can use the LO impact factor, assuming that the main energy-dependent 
NLO corrections come from the NLO BFKL sub-process rather than from the photon impact factors [26]. In order 
to simulate the NLO effects, we have used a value LOp = 0.37, shown in the dot-dashed curve of Fig. 3. The effective 
behavior obtained is given by W^^, with A = 0.29, which value is somewhat close to A ~ 0.23 — 0.29 obtained in [8]. 

It would be timely compare our estimates with the more refined works in [8] and in Ref. [28]. In the former, the most 
part of the NLO c;orrec;tion comes from considering a kinematical constraint. That is, the effects from restricting the 
range for the term of real emission in the BFKL equation, imposed upon the available phase space, corresponding to 
the requirement that the virtuality of exchanged gluons in the BFKL ladder is dominated by the transverse momentum 
squared. Such a procedure covers up to 70 % of the NLO corrections to the Pomeron intercept. As discussed above, 
our results are similar. In Ref. [28], the dipole-dipole approach is considered and predictions to double meson (light 
and heavy) production are performed. There, the result for the total cross section 77 —> J/^ J/^ is one or two order 
of magnitude below the presented here and in [8] . 

Having determined the 77 cross section we can estimate the double J/^t cross section in ultraperipheral heavy ion 
collision [Eq. (3)]. In Fig. 4 we present the energy dependence of the cross section, considering the Pb + Pb collisions 



ln"(fc7fe'^) 



2aln(W2/s) 



(16) 
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and distinct approximations for the 77 subprocess. We have found that the LLA BFKL sohition predicts large values 
of cross section for LHC energies in comparison with the Born approximation. Basically, we have that the ratio 
between the BFKL prediction and the two-gluon cross section is around 40 at ^/s = 5500 GeV. If the non-leading 
corrections are considered, as discussed above, we have that this ratio is reduced to approximately 1.7. Therefore, our 
results indicate that if the NLO corrections to the BFKL approach are account in the approximation considered here, 
a future experimental analyzes of double J/^* production in ultraperipheral heavy ion collisions could not constrain 
the QCD dynamics, since the BFKL(NLO) result is similar to the Born prediction. Only accurate measurements 
would allow discriminate between the two cases. Moreover, it is important to salient that the NLO corrections can 
be larger, implying a full suppression of energy enhancement associated with the iteration of gluons in the t-channel 
present in the QCD dynamics at high energies. However, if the energy dependence of the 77 cross section was driven 
for a large intercept, closer to the LO prediction, the analyzes of ultraperipheral heavy ion collisions can be useful. 

Using the results for the cross section shown in Fig. 4 we can estimate the expected number of events for the 
LHC luminosity. For PbPb collisions with energies of center of mass equal to = bbQQ AGeV , luminosities of 
Caa = 4.2 X lO^^cTO^^ s^^ arc used. Consequently, during a standard 10^ s/ month heavy ion run at the LHC, we 
predict approximately 30, 50 and 1100 events for Born, BFKL (MOD) and BFKL (LO), respectively. It is interesting 
to compare these predictions with the results for double J/\l/ production at LEP2, since the 77 center of mass energies 
are similar. Considering an integrated luminosity of 500p&~^ in three years and ^/s = 175GeV, almost 70 events are 
expected taking into account the non- leading corrections to the BFKL approach [8]. Therefore, we predict a large 
number of events in ultraperipheral heavy ion collisions, allowing future experimental analyzes, even if the acceptance 
for the J/^f detection being low. It should be stressed that the LHC probably will operate in its heavy ion mode only 
four weeks per year. 

Some comments related to background processes are in order here. We did not consider the additional contribution 

of the pomeron-pomeron process in the present calculations. It has been verified that such reactions would be non- 
negligible for light ions, while they are significantly suppressed for heavy ions [18]. In the particular case of the double 
heavy meson production this contribution deserves more detailed studies, since the current treatments rely on the 
Regge formalism instead of a QCD approach. An important background for the photon-photon processes are the 
photonuclear interactions, since the reactions have similar kinematics. For the process consider here, the diffr active 
J/^ production in photon-pomeron interactions [14,29] should contribute significantly. In particular, because the 
cross section for this process is large, the probability of having double (independent) production of J/\E' in a single 
nucleus-nucleus collision, associated to multiple interactions, is non- negligible [29]. However, in principle, an analyzes 
of the impact parameter dependence should allow to separate between the two classes of reactions, since two-photon 
interactions can occur at a significant distance from both nuclei, while a photonuclear interaction must occur inside or 
very near a nucleus. We salient that the experimental separation between the two classes of processes is an important 
point, which deserves more studies. An additional contribution in two photon ultraperipheral collisions is the meson 
production accompanied by mutual Coulomb dissociation. It has been estimated that these reactions increases the 
total cross section for an amount of ~ 10% at RHIC/LHC energies [30]. We disregarded this contribution in our 
calculations, since we believe that it is smaller than the theoretical uncertainty, associated, for example, to the choice 
of the running coupling constant, charm mass, etc. 

IV. SUMMARY 

Dedicated measurements have been proposed for detecting and studying the large 1/x logarithm ressumation effects 
in QCD. Experimentally establishing the BFKL effect in data is very important for the understanding of the high 

energy limit in QCD scattering. In the last years, many studies with this objective have been realized, but the 
conclusions are not unambiguous. This situation should be improved in the future with the next generation of linear 
colliders. As an alternative, in this paper we propose to investigate QCD Pomeron effects in a different context, 
namely, in photon-photon scattering at ultraperipheral heavy ion collisions. We have analyzed the potentiality of 
this process to constraint the QCD dynamics at high values of energy and provide reliable estimates for the cross 
sections concerning that reaction. In particular, we analyze in detail the double J/* production, considering distinct 
approximations for the QCD dynamics. Wc have that the LO BFKL cross section is much larger than the two-gluon 
cross section. Unfortunately, the higher order corrections of the BFKL equation are large, reducing the ratio between 
BFKL and two-gluon cross section. It makes difficult to discriminate between non-leading effects and the Born 
approximation. However, the number of events predicted for LHC is large, allowing a future experimental analyzes 
of this process as well as more detailed studies about the QCD dynamics. 

As a last comment, it is important to salient that there are limitations on the perturbative treatment that are 
intrinsic to the BFKL equation, which come from the region of very high s. The BFKL equation is known to give rise 
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to violation of the unitarity bound at asymptotically large energies. Consequently, the growth of the cross section 
predicted by the BFKL equation cannot continue indefinitely, and unitarity corrections (saturation effects) must arise 
to slow it down. They correspond to multiple bare pomeron exchanges and multi-pomeron interactions, taming the 
steep increase at high energies. We expect that the inclusion of these effects in our calculations should imply in a 
similar behavior of the non-leading BFKL predictions presented here. More detailed studies related to the unitarity 
corrections and BFKL dynamics will be analyzed in the sequel of this paper [3f], where we calculate the heavy quark 
production in ultraperipheral heavy ion collisions and investigate the saturation effects in the cross section. 
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4.21 pb 
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HERA 
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as = 0.200 


HERA 


0.97 pb 



TABLE I. The estimates for atot{^^ — > J/^ J/^) at the Born level for different Ug values. 




FIG. 1. The QCD pomeron exchange mechanism in ultraperipheral heavy ion collisions. 
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FIG. 3. The estimates for the total 77 cross section using LLA BFKL solution (solid curve) and non-leading effects 
(dot-dashed curve). 
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FIG. 4. The estimates for the double J /'i! production in ultraperipheral heavy ion collisions considering the LLA BFKL 
solution (solid curve), non- leading effects (dashed curve) and the Born approximation (dot-dashed curve). 
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